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^ . Abstract 

Consider a convex set S = {x 6 V : G(x) >r 0} where G(x) is an m x m symmet- 
C*") \ ric matrix whose every entry is a polynomial or rational function, T> C R" is a domain 

where G(x) is dchncd, and G(x) >; means G(x) is positive semidefinite. The set S is 
r \'j called semidefinite programming (SDP) representable or just semidefinite representable 

Oif it equals the projection of a higher dimensional set which is defined by a linear matrix 
inequality (LMI). This paper studies sufficient conditions guaranteeing semidefinite rep- 
el ■ resentability of S. We prove that 5* is semidefinite representable in the following cases: 
(i) T> = R", G(x) is a matrix polynomial and matrix sos-concave; (ii) T> is compact 
convex, G(x) is a matrix polynomial and strictly matrix concave on T>\ (iii) G{x) is a 
matrix rational function and q-module matrix concave on T>. Explicit constructions of 
SDP representations are given. Some examples are illustrated. 

Keywords convex sets, linear matrix inequality, matrix concavity, polynomial matrix in- 
\^ | equality, rational function, semidefinite programming, sum of squares 
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Q\ \ 1 Introduction 

p 

Suppose S is a convex set in M n given as 

S = {x G V : G(x) h 0}. (1.1) 



X 



Here T> C W 1 is a domain, and G{x) is a m x m symmetric matrix polynomial, that is, every 
entry of G(x) is a polynomial in x. The notation A y (resp. A y 0) means the matrix A 
is positive semidefinite (resp. definite). Suppose G(x) has total degree 2d and 

G{x)= ]T G a x? (1.2) 

a£N n : oiH hOn<2d 

The G a are constant symmetric matrices. The G(x) y is called a polynomial matrix 
inequality (PMI). When G(x) is linear, optimizing a linear functional over S becomes a 
standard semidefinite programming (SDP) problem. SDP is a very nice convex optimization, 
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has many attractive properties, and can be solved efficiently by numerical methods. We refer 
to [13 Ell EH [27] . There are standard packages like SDPA @], SDPT3 [25] and SeDuMi [23] 
for solving SDP problems. It would be a big advantage if an optimization problem can be 
formulated in SDP form. So we are very interested in knowing when and how the set S can 
be represented by SDP. 

One elementary approach for this representation problem is to find symmetric matrices 
A±, . . . , A n such that 



If such Ai exist, we say S has a linear matrix inequality (LMI) representation and S is LMI 
representable. Unfortunately, not every convex set in W 1 is LMI representable. For instance, 
the convex set {x 6 ]R 2 : 1 — x\ — x\ > 0} is not LMI representable, as proved by Helton and 
Vinnikov [9]. Therefore, we are more interested in finding a lifted LMI representation, that 
is, in addition to Ai, finding symmetric matrices B%, . . . ,-B/v such that 



If such matrices Ai and Bj exist, we say S is semidefinite programming (SDP) representable 
or just semidefinite representable, and (|1.3p is a lifted LMI or semidefinite representation for 
S. The variables yj are called lifting variables. Nesterov and Nemirovski [17] . Ben-Tal and 
Nemirovski [2], and Nemirovski [18] gave collections of convex sets that are SDP representable. 
Obviously, to have a lifted LMI, a convex set must be convex and semialgebraic, i.e., it can 
be defined by a boolean combination of scalar polynomial inequalities. However, it is unclear 
whether every convex semialgebraic set has a lifted LMI or not. 

When G{x) is diagonal, i.e., S is defined by scalar polynomial inequalities, there is some 
work on the semidefinite representability of S. Parrilo [21] constructed lifted LMIs for planar 
convex sets whose boundaries are rational planar curves of genus zero. Lasserre constructed 
lifted LMIs for convex semialgebraic sets satisfying certain conditions like bounded degree rep- 
resentation (BDR) [121113]. Their constructions use moments and sum of squares techniques. 
In [6], Helton and Nie proved sufficient conditions like sos-convexity and strict convexity, 
which validate lifted LMIs from moment type constructions. Later in [TJ, they further proved 
every compact convex semialgebraic set is SDP representable if its boundary is nonsingular 
and positively curved. Recent work in this area can be found in [H El El QUI Ell EH E5J ES [20] . 

We might consider to apply the existing results for the case of scalar polynomial inequal- 
ities like in [6l [7j [T2"l [T3l EH] to the case of matrix polynomial inequalities. Note 



Here polynomials pi(x) are principal minors of matrix G(x) with row (or column) index /. 
So all the results in [BJ [7J E2 E3] can be applied to study the semidefinite representability 
of S by investigating the properties of principal minors gi(x). If every gi(x) is sos-concave, 
or S is compact convex and its boundary is nonsingular and positively curved, then it is 
semidefinite representable as shown in [6l [TJ. However, these conditions are based on the 
principal minors of G(x) or the geometry of S, and hence are difficult to be applicable in 



S = {x £ R n : A) + Aixi H h A n x n >z 0}. 




(1.3) 
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many situations, especially when the matrix G(x) has big dimensions. The conditions directly 
on G(x) are preferable in practical applications. The motivation of this paper is to construct 
SDP representations for S and prove sufficient conditions directly on G(x) justifying them. 

In some particular applications, G(x) might be a matrix rational function, i.e., every 
entry of G{x) is a rational function. This is often the case in control theory. When G{x) is a 
scalar rational function, the author in [19] studied SDP representability of S. In [19], explicit 
constructions of lifted LMIs are given, and sufficient conditions validating them are proved. 
In this paper, we will construct lifted LMIs for the more general case that G{x) is a matrix 
rational function, and prove sufficient conditions justifying them. 

This paper is organized as follows. Section [2] discusses the semidefinite representation 
of S when T> = W 1 , and G(x) is a matrix polynomial and matrix sos-concave. Section [3] 
discusses the semidefinite representation of S when T> is a compact convex domain, and G(x) 
is a matrix polynomial and strictly matrix concave on T>. The case that G{x) is a matrix 
rational function and is q-module matrix concave over T> will be discussed in Section H] 

Notations. The symbol N (resp., R) denotes the set of nonnegative integers (resp., real 
numbers). For any t £ M, \t~\ denotes the smallest integer not smaller than t. The 
denotes the nonnegative orthant. For x £ M 71 , Xi denotes the i-th component of x, that is, 
x = (xi, . . . , x n ). For a £ N n , denote \a\ = a± + • ■ • + a n . For x £ W 1 and a £ N™, x a denotes 
x" 1 ■■■x"". For a, (3 £ N n , denote a < (3 if every ai < (3{. The symbol N<fc denotes the 
multi-index set {a £ N n : |a| < k}. For every integer i > 0, denotes the i-th. standard unit 
vector. The [x]d denotes the vector of all monomials having degrees at most d and ordered 
gradedly alphabetically, that is, 

[x]J = [l Xl ••• X n x\ X\X2 ■■■X 2 n ■■■ xf xf~ l X 2 •■• X n ]. 

A polynomial p(x) is said to be a sum of squares (sos) if there exist finitely many polynomials 
qi(x) such that p(x) = ^ qi(x) 2 . A matrix polynomial H(x) is called sos if there is a matrix 
polynomial F(x) such that H(x) = F{x) T F{x). For a set S, int(S) denotes its interior, 
and dS denotes its boundary. For u £ R , ||n||2 denotes the standard Euclidean norm. 
For a matrix X, X T denotes its transpose, ||^||f denotes the Frobinus norm of X, i.e., 
— \/Trace(X T X) , and ||X||2 denotes the standard operator 2-norm of X. The symbol 
• denotes the standard Frobinus inner product of matrix spaces, and In denotes the N x N 
identity matrix. For a function f(x), Z(f) = {x £ W 1 : f(x) = 0}, S/ X f(x) denotes its 
gradient with respect to x, and V xx f(x) denotes its Hessian with respect to x. 

2 Matrix sos-concavity 

In this section, assume the domain T> = M n is the whole space and G(x) is an m x m 
symmetric matrix polynomial of degree 2d. We will first construct an SDP relaxation for S 
using moments, and then prove it is a lifted LMI when G(x) satisfies certain conditions. 

A natural SDP relaxation of S can be obtained through moments. Define linear matrix 
pencils G{y) and Ad{y) as 

Q{y)= J2 y<* G -> Mv)= E 
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where G a are from (H2D and A£ } are such that 

«eN< 2d 

Since S = £ 1" : G(x) ^ 0, ^ 0} , we know 

5 = {(y ei ,...,y en ) eR n : 3xeK n ,y a = x a Vae N< 2d , 0(y) £ 0, ^(y) >= 0} . 

Here each e, denote the i-th standard unit vector whose only nonzero entry is one at index 
i. If the condition y a = x a is removed in the above, then S is a subset of 

L = {x£R n : 3y£ rCmO , y = 1, a* = y ei , • • • , x n = y Cn , Q{y) h 0, A d (y) h o} . (2.2) 

So 5 C L. Does S = L? What conditions make S = L? This section will look for sufficient 
conditions ensuring S = L. 

The matrix- valued function G(x) is called matrix concave if for all u,v £ M n and < # < 1 

G(#u + (1 - 0» ^ 0G(u) + (1 - 0)G(«). 

The matrix concavity of G(x) is equivalent to 

-V xx (£ T G(x)£) h0 V£ € K"\ Vx G R". 

Generally it is difficult to check the matrix concavity of matrix polynomials. Even for the 
simple case that they are quadratic, the problem is already NP-hard, as shown below. 

Proposition 2.1. It is NP-hard to check the matrix concavity of quadratic matrix polyno- 
mials. 

Proof. Let fn = ^m(m+l). For any symmetric matrices A±, . . . , A m 6 W nxm and B\, . . . , B m £ 
]R nxn , define matrix polynomial 

G(z) = --£(s T B i x)i4 i . 

i=i 

Then we have 

fn 

-V xx (£ T G(x)0 = J2(Z T M)Bi. 

i=l 

So G{x) is matrix concave if and only if the following bi-quadratic form 

Tn 

^2(fAiO(x T B iX ) 

i=l 

is always nonnegative. It has been proven in [16] that it is NP-hard to check the nonnegativity 
of bi-quadratic forms. Therefore, it must also be NP-hard to check the matrix concavity of 
quadratic G{x). □ 
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A stronger but easier checkable condition than matrix concavity is the so called matrix 
sos-concavity. We say G{x) is matrix sos-concave if for every £ £ M. m , there exists a matrix 
polynomial F^(x) in x such that 

- V^(£ T G(x)£) = F^xfF^x). (2.3) 

The above F^(x) has n columns but its number of rows might be different from n, and its 
coefficients of x a depend on £. 

Theorem 2.2. Suppose G(x) y for some x. If G(x) is matrix sos-concave, then S = L. 

Proof. We have already seen S C L, so it suffices to prove the reverse containment. Otherwise 
suppose L ^ S, then there must exist a point x 6 L/S*. Since S is closed and convex, by 
Hahn-Banach Theorem, there exists a supporting hyperplane H = {i€ M n : a T x > b} D 5 
such that a T u = 6 for some u £ S and a r i; < 6. Consider the linear optimization 

min a T x subject to G(x) y 0. (2.4) 



Clearly u is a minimizer and b is the optimal value. The optimization (|2.4p is convex. The 
existence of x with G(x) y 0, i.e., the Slater's condition holds, implies there exists a matrix 
Lagrange multiplier A y such that 

A.G(n)=0, a = V x {A»G(x)) 

The value and gradient of a T x — A • G(x) — 6 vanish at u. Then, by the Taylor expansion at 
u, we have 

a T x - A • G{x) - b = (x - uf (J J -V M (A • G(u + s(x - u)) ds dt^J (x - u). 

Since A y 0, there exist vectors such that A = X)fc=i ^ k \\^) T . So we have 

a T x - A • G(x) - b = 



fc=i v - / ° - 70 



- V ra ((A (fc) ) T G(u + s(x - n))A {fc) (x — u . 



Since G(x) is matrix sos-concave, by Lemma 7 in [6], we know each summand in the above 
must be sos. Thus a T x — A • G(x) — b must also be an sos polynomial of degree 2d. So there 
exists a symmetric matrix W y such that the identity 

a T x — A • G(x) - b = [x] T d W[x] d = W . {[x] d [x]J) 
holds. By definition of matrices in (|2.ip . we have 



a T x — A • G(z) - 6 = W • [ ) . 

v aeN< 2d 
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Figure 1: The drawn body is the convex set in Example 12.31 



Since x 6 L, there exists y such that x = (y ei , ■ ■ ■ ,Ve n ), Q{y) ^ and Ad{y) h 0. So in the 
above identity, if we replace every monomial x a by y a , then 



a x — A • Q(y) - b 



or equivalently 



a T x — b 



w* [J2 =w.A d (y), 

\\a\<2d J 

= AmG(y) + W*A d (y). 



Since A, 0(y), W, Ad(y) t 0, we must have a T x — b > 0, which contradicts the previous 
assertion that a T x — b < 0. So S = L. 

Example 2.3. Consider the set {i e I 3 : G{x) >z 0} where 



□ 



G(x) 



-x{- 2xi 

1 + X\X2 
XlX 3 



1 + XiX 2 



2x1 



1 + x 2 x 3 



1 + x 2 x 3 

2 ^"-^2 



The plot of this set is in Figure [TJ The Hessian — V X x(£, T G(x)£,) is positive semidefinite for 
all £ € M 3 . This is because 



-V ra (£ T G(x)0 



-66 
tl + 26 2 

3 



-66 -66 6 2 + 2£ 2 



-66 
-66 
,3 + 26, 



h V£ G 



which is due the fact that the bi-quadratic form x T ( — V xx (^ T G(x)^)^x in 

afci + + 4tl + 2 (x\& + Ail + z§£?) " 2(xix 2 66 + ^1^366 + x 2 x 3 66) 
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is always nonnegative, as shown by Choi [3]. So this set is convex, and by Theorem 12.21 a 
lifted LMI for it is 



x£l 3 : 3y a (a £ N< 2 ) such that 
2 - 2/200 - 2y o2 1 + yno 2/101 
1 + 2/no 2 - 2/020 - 22/200 1 + 2/on 
2/ioi 1 + 2/on 2 - 2/002 - 22/020 



>=0, 



1 


xx 


»2 


ai3 


ail 


2/200 


2/no 


2/ioi 


ai2 


2/no 


2/020 


2/on 


ai3 


2/ioi 


2/on 


2/002 



>- 



There are in total 6 lifting variables. The plot of the x in the above coincides with the convex 
body of Figure [H which confirms this lifted LMI is valid. □ 

The matrix sos-concavity condition requires us to check the Hessian 

is sos for every £ £ M m . This is almost impossible since there are uncountably many cases. 
However, a stronger condition called uniformly matrix sos-concave is 



-V xx .(£ r G(x)0 = FfoxY F(£,x), 

where F(£,x) is now a matrix polynomial in joint variables It is easier to check. The 

uniformly matrix sos-concavity can be verified by solving a single semidefinite programming 
feasibility problem. Clearly the following is a consequence of Theorem | 



T ; 



Corollary 2.4. Suppose G(x) >- for some x . If G{x) is uniformly matrix sos-concave, 
then S = L. 

It should be pointed out that when G(x) is matrix sos-concave, it is not necessarily that 
G(x) is uniformly matrix sos-concave. For a counterexample, consider the G(x) defined in 
Example 12.31 the quadratic matrix polynomial —V X x{£ T G(x)£) there is not sos in £, as shown 
by Choi [3]. Now let us see an example of uniformly matrix sos-concave G(x). 



Example 2.5. Consider the set {i £ I 2 : G{x) y 0} where 



G(x) 



2x\ 



x\x2 



kx\x\ 



X1X2 



1x^ 



3 — x\x2 



- XIX2 

2J2 



\X\X 1 2 



A plot of this set is in Figure [2j The above G(x) is uniformly matrix sos-concave because 

-y xx (C T G(x)0 =H 1 + H 2 + H 3 + H 4 , 

1 T 



Hi 



26 x\ + &X2 
H, = 



26^1 +C2X2 
1£\x 2 + £2x1 

£1X1 £2X2 £2X1 
£2x1 £1x2 £2x2 



H2 = 8fe 2 + £2 

i\x\ &X2 &xi 
£2x1 £1x2 £2x2 



011x2 



X\X 2 

™2 
,L ■-) 



-1 T 



H4 = 2(((fxf + f2xl) 



1 
1 



+£l 



2x\ + \x\ 




7 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 

-0.6 -0.4 -0.2 0.2 0.4 0.6 

Figure 2: The shaded area is the convex set in Example 12.51 and the thick curve is its 
boundary. 

So this set is convex, and by Corollary 12.41 a lifted LMI for it is 

i£l 2 : 3 yij (0 <i,j < 4) such that 

\2 - 2(y 40 + 2y 22 + 2/04.) 3 - (y sl + y 13 ) 
3 - (2/31 + yrs) 5 - (2/40 + 32/22 + Voa) 

There are in total 12 lifting variables yij. The plot of x in the above coincides with the shaded 
area of Figure [21 which confirms this lifted LMI is valid. □ 

3 Strict matrix concavity 

In this section, assume S = {x £ T> : G(x) >z 0} where 

V = {xeR n : 9l (x) >0,...,g m (x) > 0} 

is a domain defined by multivariate polynomials gi(x), . . . , g m (x). When T> is compact convex 
and G(x) is strictly matrix concave on V, we will show S is semidefinite representable, and 
a lifted LMI for it can be explicitly constructed. 

Like in the previous section, a natural SDP relaxation of 5 can be constructed by using 
moments. Let go(x) = 1, and set d = max {deg(G(x))/2, [~deg(<?fc)/2] , k = 0,1, . . . ,m} . For 
every integer N > d and k = 0, . . . ,m, define symmetric matrices such that 

9k{x)[x] N ^ dk [x] T N _ dk = x ^ B k!p- ( 3 - 1 ) 

/3eN< 2 ]v 




1 X\ x 2 

xi y2o 2/n 

X2 2/11 2/02 

2/20 2/30 2/21 

2/11 2/21 2/12 

.2/02 2/12 2/03 



2/20 2/ii 2/02 

2/30 2/21 2/12 

2/21 2/12 2/03 

2/40 2/31 2/22 

2/31 2/22 2/13 

2/22 2/13 2/04 



s 



This determines B^O uniquely. Then define linear matrix pencils B^ (y) as 

B ( k N) (y)= E yf> B k$> fc = o,i,...,m. 

/3GN< 2Ar 

Clearly 5 can be equivalently described as 

q \ i , cm n ^xeR n ,y a = x a VaeN< 2ff) 1 

5= l (yei '"" yejGM : g(y)±0,Bl N \y)yO,k = 0,...,m]- 

If the condition y a = x a is removed in the above, then S is contained in the set 
t f rD » q ^m>( n+2N ) yo = l,xi = y ei ,...,x n = y en , \ 

iw = |xel :3!/ €R( . J, 5(!) , 0iS l») (j)Ht = 0i . ..,„)■ < 3 ' 2 > 

So we have 5 C L^r for every N > d. Clearly L/v+i C L^r, because -Lat+i is a restriction of 
Ltv. Thus it holds the nesting containment relation: 

La 2 ■ ■ ■ 2 -^n 5 £jv+i 5 • • • 5 <S. 

Does there exist a finite iV such that S = L^l What conditions on S make it true? This is 
the main topic of this section. 

Semidefinite representation of S is closely related to linear functionals nonnegative on S. 
For a given ^ £ G R n , consider the linear optimization 

min £ T x subject to G(x) >z 0. (3-3) 

x&T> 

It always has a minimizer u G dS when S is compact. If S C int(T>), then u G int(T>), and 
hence the constraint 1 £ P in (|3.3p is not active. If further there exists x E T> such that 
>- 0, i.e., the Slater's condition holds, and G(x) is matrix concave in T>, then there 
exists a matrix Lagrange multiplier 0^A6 M raxm such that 

A.G(n)=0, e = V x (A.G(x))\ x=u . (3.4) 
Thus, by its Taylor expansion of at u, we know i T (x — u) — A • G{x) equals 



(x — u) 



J J -V xx (A»G(u + s(x-u)))dsdtj -(x-u). 



(3.5) 



H(u,x) 



If the above matrix polynomial H(u, x) has a weighted sos representation in terms of G(x) 
and gi(x), then we can also get a similar one for l T [x — u). For this purpose, we need some 
assumptions on T> and G(x). 

Assumption 3.1. G(x) is matrix concave on T>, and G(x) satisfies 

-^V„(A.G(ti))^0 VuG dS, VO/A^O. 
A F 
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Theorem 3.2. Suppose 5 C int(T>), T> is compact convex, and there exists x G T> such that 
G(x) y 0. If Assumption \3. 1\ holds, then S = Ln for all N big enough. 

Proof. For a matrix polynomial G(x) given in (jl,2p . we define its norm ||G(x)|| as 



\\G(x)\\ = max —. ' \\G, 



aGN< 2d \Ot\\ 

Since T> is compact, there must exist A > such that 
1 



a 2- 



|tf(u,x)||<A VAGl nxn , V« G 95, ViGD. 

Here H(u,x) is defined in (13. 5p . Assumption 13.11 implies H(u,x) >- for all u € (95 and 
x € T>. This is because otherwise if H(u, x) is not positive definite, we can find / u 6 R n 
such that v T H(u,x)v = 0, i.e., 

J J v T ( - (A • GO + s(x - u)) ds dt = 0. 

Since G(x) is matrix concave on the convex domain T>, we must have 

V xx (A*G(u + s(x-u)))v = V0<s<i<l. 



T 



In particular, we get v T y — V X x (A • G{uj) J v = 0, which contradicts Assumption 13. 1[ There- 
fore, by the compactness of dS and V, there exists 5 > such that 

-H(u,x)hSI n VuedS, Vxe£>, VO^A^O. 



||A||. 

By Theorem 27 in [BJ, there exists an integer iV* such that for every ^ A ^ and it € OS 1 , 
there exist sos matrices Fq(x),Fi(x), . . . , F m (x) satisfying 

j m 
|-^#(u,x) = ^ 5fc (x)F fe (x), (3.6) 

deg{F k ) + 2d k <2N*, k = 0,,...,m. 

Now we claim S = Ln*. Since 5 C Ljv», we need to show Ltv* C S. Suppose otherwise 
there exists x € Ln*/S. Since T> is compact convex and G(x) is matrix concave on T>, 5 is 
closed and convex. By Hahn-Banach Theorem, there exist ^ I E M n and u € dS satisfying 

f{x-u)>0 VieS, ^ T (x-u)<0. 

Consider the linear optimization (|3.3|) with this The point u € 55 is a minimizer. Since 
5 C int(T>), the constraint x £ 2? is not active in (|3.3p . The existence of x € T> with G(x) >- 
(the Slater's condition holds) implies there exists A >z satisfying (13. 4ft . From (13. 5p and ()3.6p . 
we know there are SOS polynomials po(^)i Pi(^) 5 • • • , Pm(x) such that 

m 

f(x-u) = A.G(x) +^2p k (x)g k {x), 

k=0 
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and they satisfy 

deg(p fc ) + 24 < 2iV*, fc = 0,l,...,m. 
So there are symmetric matrices Wq, W±, . . . , W m y such that 



f(x - u) = A • G(x) +^5fe(^)N^-d fe ^[ a; ]Ar*-d fc - 



fc=0 



By definition of matrices B^ in (13. 1|) , it holds the identity 

m / 

By the choice of x, there exists y such that x = (y ei , . . . , y e „), £?(y) ^ 0, and every Bf? (y) y 0. 
So in the above, if every monomial x a is replaced by y a , then 



^ T (£ - u) = A . + J2 w k • ^f* (y) > o, 



fc=0 

which contradicts £ T (x — u) < 0. So we must have S = Ln*. 

For every N > N* , the relation S Q Ljy Q Ln* implies S = Ljy. So the theorem is 
proven. □ 

Assumption 13.11 requires us to check — V xa; (A • G(u)) >- for every nonzero A y and 
u € dS, which is sometimes very inconvenient. However, Assumption 13.11 is true if G(x) is 
strictly matrix concave on T>, that is, for every ^ £ £ M m , the Hessian — V a;a: (^ T G(a;)^) >- 
for all x € D. So the following is a consequence of Theorem 13.21 



Corollary 3.3. Suppose S C int(T>), T> is compact convex, and there exists x € V such that 
G(x) y 0. If G(x) is strictly matrix concave on T>, then S = for all N big enough. 



4 Rational matrix inequality 

In this section, assume S = {x € T> : G(x) y 0} is defined by a matrix rational function G(x), 
i.e., every entry of G(x) is a rational function. Suppose G(x) is matrix-concave on T>. As 
before, the domain T> = {x € M. n : g\{x) > 0, . . . ,g m (x) > 0} is still defined by multivariate 
polynomials. When G{x) is a scalar rational function, the author in [19] studied the SDP 
representation of S. This section discusses the more general case of G{x) being a matrix. 
We will first construct an SDP relaxation for S, and then prove it represents S when G(x) 
satisfies a certain condition. 

Suppose the matrix rational function G{x) is given as 

v ' aSN": |a|<deg(G) 

where F a G R mxm are symmetric matrices, den(G) is the common denominator of G(x), and 
deg(G) is the degree of G(x). Let p(x),q(x) be two given polynomials which are positive in 



11 



int(V). We say G(x) is q-module matrix concave over T> with respect to (p, q) if for every 
£ € M m , there exist SOS polynomials ai t j(x,u) such that 



• U T G(u)Z + (V x eG(x)0 T (x-u)- eG(x)t) 

\ , x=u , / 



E 0*0*0 Y.9j(u)<7ij(x,u) 

*=0 \jf=0 



(4.2) 



is an identity in (a;, u). In the above goC^) = 1- If (PiQ) can De chosen as (den(G),den(G) 2 ), 
then we just say G(x) is q-module matrix concave over T>. The condition (|4.2p is based on 
Putinar's Positivstellensatz [22]. Clearly, if G(x) is q-module matrix concave over T> with 
respect to (p,q), then it must also be matrix concave over T>. 

Now we turn to the construction of lifted LMI for S. Assume G(x) is q-module matrix 
concave over T> with respect to (p, q) and (14. 2D holds. Set integer 



d = max <^ max \\ deg^^cJi,)] , \ deg(G) \ . 

\0<i,j<m ) 



(4.3) 



For i = 0, . . . , m, define matrices P<^ , Qa^ such that 

f^T Md-<fi Nd-dj = E Q^X a + E Pa^^by- (4,4) 

a&N n :\a\ + \LE(p)\<2d l3eN< 2d -.p<LE(p) X ' y J 

Here LE{p) denotes the exponent of the leading monomial of p(x) in the alphabetical ordering 
(xi > X2 > ■ • • > x n ). Let y be a vector indexed by a S N n with \a\ + \LE{p)\ < 2d, and z 
be a vector indexed by (3 € N<2d with /? < LE(p). Define linear matrix pencils 

Qi{y,z)= E Qa ] y a + E p « 39 > i = o,i,...,m. (4 5) 

a£N n :\a\ + \LE{p)\<2d f3m< 2 d:/3<LE(p) V ; 



Here P%> ,Ql> are defined in (|4.4j) . Suppose G(x) is given as 

,/3 



,0 n(*) 

C 

g(x)= e F « )xa+ £ F f ] 

«GN n :|a|- 

then define linear matrix pencil 



p(x) 



aGN":|a| + |L£;(p)|<2d /3GN< 2d :/3<L£:(p) 

Clearly, 5 can be equivalently described as 

[ h 0, Qi(y,z) h 0, * = 0, ...,m J 

If we remove y Q = x a , zp = in the above, then S is a subset of 

r -J-cion. 3y,z, y = = 7/ ei ,...,x n = y en , 1 , . 

V_ l e • F(y,z)>rO J Q i (y,z)>:0,t = 0,... J m J" ^ 

So S C L qm0 (i. We are interested in conditions making 5 = L qmof i. 
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Lemma 4.1. Assume S C int(T>) and there exists x G T> such that G(x) y 0. Suppose the 
matrix rational function G(x) is q-module matrix concave over T> with respect to (p,q). If 
v G dS, q(v) > 0, den(G)(v) > 0, and a T {x — v) > for all x G S, then 

m 

p(x) ■ (a T (x - v) - A • G(x)) = y^gi{x)ai{x) 

i=0 

/or some symmetric matrix A ^ and sos polynomials o~j(x) mi/i deg(a;,o"i) < 2d. 
Proof. Consider the linear optimization 

min a T x subject to G(x) y 0. 

x(zD 

The point u € OS is an optimizer. Since den(G)(v) > 0, G(x) is differentiable at v. Since 
S C int(D), the constraint x G 2? is not active at v. Because of the existence of x G £> with 
G(x) y (the Slater's condition holds) and G(x) being matrix concave on T>, there exists a 
Lagrange multiplier A y such that 

a = V x (AmG(v)), A«G(u) = 0. 

Hence, we get the identity 

a T (x - v) - A • G(x) = A • G(v) + V X (A • G(v)) T (x - u) - A • G(x). 
Since A y 0, we have decomposition A = X^i=i A^(A^) T . Then it holds 

a T (x - v) — A • G(x) = 

{(X {l) ) T G(v)X^+V x ((X^) T G(v)X^f(x -v)- (A«) T G(x)A«}. 

So the lemma readily follows the q-module matrix concavity of G{x). □ 

For a function f(x), denote by 2(f) its real zero set, i.e., Z(f) = {x G M n : /(x) = 0}. 

Theorem 4.2. Assume S is closed and convex, S C int(T>), G(x) >- /or some ieP, and 

dim(Z( det(G)) n Z{den(G)) HdS) <n-l, dim (2(g) n dS) < n - 1. 

If G(x) is q-module matrix concave over T> with respect to (p,q), then S = L qm od- 

Proof. Since S C L qmo ^ it is sufficient for us to prove the reverse containment. For a 
contradiction, suppose there exists x 6 L grnoc i/S and (j/, 2) such that 

x = (yei,---,ye n ), F(y,z)h0, Qi(y,z)y0,i = 0,...,m. 

Since S is convex and closed, by Hahn-Banach Theorem, there exists a supporting hyperplane 
{a T x = 6} of 5 such that a T x > b for all x G S and a T x < b. Let v G <9<S be a minimizer 
of a T x on S". Since dim (det(G) n Z(p) (1 dS) < n - 1 and dim (Z(q) n <9S) < n - 1, by 
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continuity, the supporting hyperplane {a T x = b} can be chosen to satisfy deg(G)(v) > 0, and 
q(y) > 0. By Lemma 14.11 we have 

rn / \ 

a T (x-v)=A»G(x) + y^f(a i (x) (4.7) 

for some sos polynomials (Ji(x) such that every deg(gj<7j) < 2d. If we write &i{x) as 

= [x]d-diWi[x]d-di 
for symmetric Wi ^ (i = 0, 1, . . . , m), then identity (|4.7|) becomes 

a T (x - v) = A . G(x) + £ f ^ Wd-* * Wi 



i=0 



= a-g(x) + e E cgw E 

t=0 \QeN™:|a|+|LE(p)|<2(i f3eN n :/3<LE(p) ^ V '/ 

In the above identity, if we replace every x a by y a and by zg, then 

a T x - b = A • F(y, i) + E *) * W i > °> 

i=0 

because all A, F(y, z), Qi(y, z),W{ y 0. This contradicts a T x < b. So S = L qmo a- D 



The condition of q-module matrix concavity requires checking (|4.2j) for every £ £ R n . In 
many situations this is almost impossible. However, if we consider £ as an indeterminant, 
then one sufficient condition guaranteeing (14. 2|) is 



(4. 



p(x)q(u) ■ U T G(u)£ + (V x eG(x)0 T (x-u)- eG(x)d) 

\ x=u / 



Y,9i(x) Y^9j(u)aij(x,u,0 , 

i=0 \j=0 J 

where every o~ij(x,u,£) is now an SOS polynomial in (x, u, £). If G(x) satisfies condition 
(|4.8|) . we say G(x) is uniformly q-module matrix concave over T> with respect to (p,q). If in 
addition p = den(G)(x), q = deg{G) 2 {u), then we just say G(x) is uniformly q-module matrix 
concave over T>. So the following corollary follows Theorem 14.21 



Corollary 4.3. Assume S is closed and convex, S C int(T>), G(x) y for some x € T>, and 

dim(Z(det(G)) n Z(den(G)) n dS) < n - 1, dim(Z(g) n 95) < n - 1. 
If G(x) is uniformly q-module matrix concave over V with respect to (p,q), then S = L qmo< i- 

Now we give some examples on how to apply Theorem 14.21 and Corollary 14.31 to get an 
SDP representation for S when G{x) is a matrix rational function. 
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, 1 1 1 1 1 1 

2 4 6 8 10 12 



Figure 3: The shaded area is the convex set in Example I4.41 and the thick curve is its 
boundary. 



Example 4.4. Consider the set {x 6 : G(x) y 0} where 
G(x) 



7 - xi + 2x 2 
5 



11 



5 


1 




- x 2 _ 


X\X 2 


xl X 2 _ 



Its domain T> is the nonnegative orthant. A plot of this set is in Figure El The above G(x) 
is uniformly q-module matrix concave over M.^_, because 

xix 2 uj u\ ■ U T G{u)Z + (V x eG(x)0 T (x-u)- eG(x)£ 

\ x—u 
x 2 ul(ui£, 2 ~ X1& + 1*1X2^1 - u 2 xi£i) 2 + xiul^(u 2 ~ X 2 ) 2 . 

Hence, this set is convex, and by Corollary 14.31 it has the following lifted LMI 
\7 



xi + 2x 2 
5 



x 2 




"0 1 

1 x l 

1 x\ x 2 

xi y 20 yu 

1 xi x 2 y 20 yn 1/02 
_0 x 2 yn y 02 



5 

11 - x 2 




ZlO + ^03 
Z02 



Z02 



h 0, 



+ 



2 210 Z01 2 20 Z Q2 



210 Z20 

z i 

^20 2 30 



2Q2 





202 




203 



230 

Z 03 

z m 



z 04 



ho, 



"0 





1" 




210 


220 


0" 




"0 


1 


0" 




201 





202 








Xl 


+ 


220 


230 





ho, 


1 


Xi 


x 2 


+ 











1 


Xl 


x 2 _ 




_ 













x 2 


0_ 




.202 





203. 



>- 0. 



The plot of x in the above coincides with the shaded area of Figure [3l which confirms this 
lifted LMI is valid. 
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Figure 4: The shaded area is the convex set in Example I4.5I and the thick curve is 
boundary. 



Example 4.5. Consider the set {x € M? : G(x) y 0} where 



G(x) 



1 — 2x1 ~ 2 x i x 2 



x\ 



x 2 



1 



T 2 
X l. 



4 


"-1 l " 


+ 2 
X-^ ~\~ X<2 


l -l 



A plot of this set is in Figured! Its domain is the whole plane R 2 , and the above G(x) 
uniformly q-module matrix concave over M 2 , because 

11*11! • INI! • U T G(u)Z + (V x eG(x)0 T (x-u)- eG(x)A = 

\ x=u / 

(g/?) ■ (6 - 6) a + • NIS • II* - u||i ■ g 

where the polynomials fi are given as below 

fl = —U\U2x\ — U\U2x\ + U1U2X2 + u\u2X\, = 775 ( _ U2 X 2 + U 2 X 2 ~ u\x\ + u\x\), 

h = ~U\U2X\ + U\U 2 x\ + U\u\x 2 - ufu2Xl, f 7 = -2u\U2X\X 2 + UlU^Xl + ufu2X 2 , 
h = -^(-U2 x l x 2 + U%X! - u\x\X2 + u\x 2 ), h = A. x \ ~ U l x 2, 

h = -^m{u\x\x 2 - u\x\ - u\x x x 2 + ufx 2 ), fg 

h = -^( U 2 X 2 - U 2 X 2 ~ u\ x \ + u l x l )> 



-U1U2X1 + u\u2X 2 - 



So this set is convex, and by Corollary 14.31 a lifted LMI for it is 



1 - 2y20 - 2yn - y 02 - z 04: y 2 o + z a 

V20 + zqa 1 - U20 ~ zqa 



y 0, 
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1 





" 




" ^00 


ZiO 


ZQ1 


~2 2 


•-'11 


202 





1 





Xl 


X2 









-ZQ2 


Zll 


-^12 


-203 


212 











XI 








+ 


zoi 


Zll 


Z02 


-203 


212 


203 


1 


X\ 


X2 


2/20 ~ 2/02 


2/n 


2/02 


-Z 2 


-Zl2 


-Z03 


204 


-213 


-2 / 










2/n 


2/02 







Zll 


~Z03 


Z\2 


-213 


-2 4 


213 











2/02 





_ 




_ ZQ2 


Z\2 


2Q3 


-204 


213 


204 



h o. 



The plot of x in the above lifted LMI coincides with the shaded area of Figure [U which 
confirms it is valid. 



5 Conclusions 

This paper gave explicit constructions of SDP representations for the set S = {x € T> : 
G(x) y 0} when G{x) is a matrix polynomial or rational function, and proved sufficient 
conditions validating them. These conditions are generally based on the matrix concavity of 
G(x). We would like to point out that G(x) might not be matrix concave when S is convex. 
For instance, for the quadratic polynomial matrix 



Q{x) 



X1X2 + 2 X\X2 

X1X2 X1X2 — 1 

Xl + X2 



the matrix inequality Q(x) y defines the convex set {x € : X1X2 > 2}, but Q(x) is not 
matrix concave on In such cases, the lifted LMIs constructed in this paper might no 
longer represent S. It is an important future work to find SDP representations for S when 
G(x) is not matrix concave. On the other hand, if S is compact convex and its boundary 
is nonsingular and positively curved, it is shown in [7] that S has a lifted LMI. But it 
is not clear how to find an efficient one. This is another interesting future work. More 
recent results on semidefinite representability of convex semialgebraic sets can be found in 

Acknowledgement The author would like to thank Bill Helton for his fruitful suggestions 
on preparing this paper. 
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